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Quantum chaotic dynamics is obtained for a tight-binding model in which the energies of the atomic 
levels at the boundary sites are chosen at random. Results for the square lattice indicate that the 
energy spectrum shows a complex behavior with regions that obey the Wigner-Dyson statistics and 
localized and quasi-ideal states distributed according to Poisson statistics. Although the averaged 
spatial extension of the eigenstates in the present model scales with the size of the system as in the 
Gaussian Orthogonal Ensemble, the fluctuations are much larger. 



Recent advances in nanotechnology have made possi- 
ble the fabrication of devices in which carriers are mainly 
scattered by the boundaries and not by impurities or de- 
fects located inside them As these devices, which 
are commonly referred to as quantum dots, resemble 
quantum billiards, the interest in the latter has increased 
considerably in the last few years. Although chaotic 
billiards have been intensively investigated in the last 
twenty years , the behavior of their quantum analogues 
has not yet been fully characterised. Some general char- 
acteristics of quantum chaotic systems have, however, a 
wide acceptance. It has been shown, for instance, that 
the quantum counterparts of billiards having chaotic tra- 
jectories, have an energy spectrum which obeys Wigner- 
Dyson statistics j§-@]. This is the case of the stadium 
and Sinai's billiards |7|-|l0|] . On the other hand, it is com- 
monly believed that there is a perfect mapping of quan- 
tum chaotic billiards into the more general and inten- 
sively investigated problem of random matrices (|,||,^],|IlJ . 

The purpose of this Letter is to present a new model 
of a quantum chaotic billiard and investigate its spectral 
statistics and its relationship with random matrices of 
the Gaussian Orthogonal Ensemble (GOE). The model 
is a practical implementation of surface roughness and 
its main characteristics are the following. The quantum 
system is described by means of a tight-binding Hamil- 
tonian with a single atomic level per lattice site in which 
the energies of the atomic levels at the surface sites S are 
chosen at random, namely, 
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where the operator c\ destroys an electron on site i, and 
Vij is the hopping integral between sites i and j (the sym- 
bol < ij > denotes that the sum is restricted to nearest 
neighbor sites) We take V. tj = V = -1. The en- 

ergies of the atomic levels at the boundary sites, e,, are 
randomly chosen between — W/2 and W/2. Calculations 
have been carried out on L x L clusters of the square 
lattice of sizes up to L — 200. Schwarz algorithm for 



symmetric band matrices |]13| was used to compute the 
whole spectrum including eigenvectors for L < 64 and 
the complete set of eigenvalues for L < 170. Instead, for 
larger matrices individual eigenvalues and eigenfunctions 
were obtained by inverse iteration . 
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FIG. 1. Variance of nearest-level spacings in the whole en- 
ergy spectrum. The circle size is proportional to the actual 
size of the system (L — 15, 30, 60, 110 and 150). The horizon- 
tal chain line indicates the value of the variance corresponding 
to the Wigner-Dyson distribution (0.286). Inset: Distribution 
of nearest level spacings in a 170 x 170 cluster for energy levels 
between -4 and -3.7 (empty circles) and -1.2 and -0.9 (filled 
circles) ; for the sake of comparison the Wigner-Dyson and the 
Poisson distributions (continuous lines) are also shown. The 
disorder parameter is W — 2. 

The basis for expecting chaotic behavior in this model 
lies on the fact that the shift of the eigenvalues pro- 
moted by the perturbation at the surface would be about 
(W/ V3)L~ 3 / 2 if first order perturbation theory were ap- 
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plicable, which is larger than the average level separation 
(~ 8£~ 2 ). This implies a mixture of ~ W\f~L ideal eigen- 
states to form a given wavefunction of the perturbed sys- 
tem. A similar reasoning suggests quantum chaos for our 
model in any dimension greater than 2. Note also that, 
in contrast with standard chaotic billiards, our model 
has two length scales, namely, the size L and the lat- 
tice constant a. Thus, even in the the macroscopic limit 
(L/a — > oo), microscopic roughness remains, and is con- 
sequently felt by quantum particles, i.e., by particles of 
wavelengths of the order of a. 




FIG. 2. Probability amplitude of the eigenstate at 
E ~ —3.3 of an 200 x 200 cluster of the square lattice with 
the energies of the atomic levels at the boundary sites ran- 
domly chosen between -1 and 1. The probability amplitude 
is roughly proportional to the darkness. 

According to the theory of random matrices and the 
numerical results for the stadium and Sinai's billiards, 
a clear hallmark of chaotic behavior is a level separa- 
tion statistics of the Wigner-Dyson type. Fig. 1 shows 
the variance of the nearest level spacings distribution for 
several cluster sizes of the billiard described by Hamil- 
tonian (1) in the full energy spectrum. Away from the 
bottom of the band the levels are distributed according 
to Wigner-Dyson statistics (this is explicitly shown in the 
inset of the figure). The behavior is, nonetheless, some- 
what different depending on the particular energy region. 
In fact, whereas for energies in the range [—2, —0.5], the 
variance is close to that of the Wigner-Dyson distribution 
(0.286) even for small clusters, away from that region the 
variance tends to 0.286 as the size is increased. On the 
other hand, near the band edges the variance of the dis- 
tribution clearly tends to 1 (uncorrelated levels) as the 
system size increases, while the distribution approaches 



Poisson distribution (see inset of Fig. 1). A similar behav- 
ior was obtained by Pavloff and Hansen in their study on 
the effects of surface roughness on metallic clusters [ p^[ : 
in the bottom of the spectrum the de Broglie wavelength 
is large, only the averaged disorder is felt and the pertur- 
bation is accordingly small. These features differ from 
those of random matrices which show a spectrum char- 
acterised by the Wigner-Dyson statistics throughout the 
whole energy range, suggesting that, at least at a meso- 
scopic level, random matrices and chaotic billiards are 
not equivalent. 
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FIG. 3. Spatial extension of all eigenstates corresponding 
to a realization of Hamiltonian (1) for a 64 x 64 cluster and 
two values of the disorder parameter, namely, W = 2 (top) 
and W = 5 (bottom). 

The spatial dependence of the probability density of an 
eigenfunction in the energy range where chaotic behavior 
is expected is illustrated in Fig. 2. The pattern shows 
speckles with valleys in between, reminiscent of the pat- 
terns found in billiards showing quantum chaotic behav- 
ior B and in prelocalized states in 2D . In comparing 
our results with those for prelocalized states it should be 
noted that in the latter case disorder is present at all 
lattice sites of the system, whereas in our model disor- 
der is restricted to the surface (it is in this sense that 
our model can be called a billiard). Another important 
difference, which is a consequence of the previous one, is 
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the fact that in the present case Anderson localization 
should not be expected. 

To further investigate the nature of the states in the 
different energy ranges we have calculated the participa- 
tion ratio, which is a good measure of the spatial exten- 
sion of a given eigenstate. Inverse participation numbers 
are defined as the moments of the distribution function 
of the local weights of the eigenstates, namely, 
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8=1 



(2) 



where a a i is the amplitude of the crth eigenstate at site 
i, i.e., \(f> a >= X)i a c"K >■ Then the participation ratio, 
P a , is given by the inverse of the second moment defined 



by (2): 



Pa — 



(3) 



P a is interpreted as the number of lattice sites covered 
by the eigenstate a. Fig. 3 shows P a in the whole energy 
spectrum for two values of the disorder parameter W . 
For low disorder (W = 2) the energy levels close to the 
band edges and part of the states close to the band center 
are quasi-ideal. The probability amplitude of eigenstates 
in that energy region is much like that of eigenfunctions 
in the fully ordered system p7[ . These states appear in 
energy regions in which the wavefunctions of the ordered 
cluster have small weight on the surface layer (at band 
edges all states have a small surface sensibility whereas 
close to the band center some of the states show small 
amplitudes at the surface) whenever the value of W is 
small enough to allow them to keep their unperturbed 
characteristics. In any case, what matters in the present 
analysis is that the size of the regions of appearance of 
quasi-ideal states diminishes both with the increase of L 
and the strength W of surface disorder. Results of Fig 
1. are consistent with this qualitative analysis. The be- 
havior is even more complex for large W (see bottom of 
Fig. 3). Bona fide exponentially localized states appear 
outside the band (see circles close to the x-axis in Fig. 3). 
Quasi-ideal states still appear at bandedges and at the 
band center but in a reduced amount. Other states show 
now a small spatial extent due to their character of bulk 
states resonating with a particular surface impurity. Ac- 
tually, we have collected a large sample of states showing 
different characteristics and allowing therefore for differ- 
ent names. To the best of our knowledge this complex 
behavior has not been pointed out in previous discussions 
of quantum chaotic billiards. In the limit of infinite disor- 
der we expect a rather simple scenario in which bulk and 
surface are decoupled, and, consequently, ordered states 
would lie on bulk sites whereas localized states would be 
located at surface sites. This is a trivial limit and the 
most interesting situations are of course expected for fi- 
nite W values. 



It is interesting to note that a calculation for random 
matrices similar to that shown in Fig. 3 gives an almost 
energy independent distribution with a finite width. An 
interesting question is how this width (or, more precisely, 
this standard deviation) evolves with the size of the sys- 
tem for both random matrices and the present billiard. 
In the case of random matrices of dimension N = L 2 
the results for the participation ratio averaged over the 
whole energy range (P =< P a >) in clusters of sizes 
L = 4, 8, ...,44, 48, can be accurately fitted by, P = 2.2 + 
0.33L 2 , and those for the relative standard deviation of 
the distribution by ajP = 1.64/ L - 3.05/L 2 (see Fig.4). 
Thus, in the assymptotic limit the ratio a /P behaves as 
1/L. On the other hand, results for the present model 
of quantum billiard obtained for L = 4, 8, ...,60,64 and 
W = 2 lead to the following fittings: P = 0.7LL + 0.33L 2 
and ajP = 0.074+1.98/L-11.8/ J L 2 +24.5/ J L 3 (see Fig.4). 
Averaging sets always include more than 8000 eigen- 
states. These results indicate that in the macroscopic 
limit a/P is a constant, suggesting that fluctuations of 
the spatial extension of wavefunctions are much larger 
than for random matrices. The reason for this significant 
difference should be the surface resonances and quasi- 
ideal states found in the present chaotic billiard which 
seem to determine the assymptotic behavior of fluctua- 
tions in this system. On the other hand we note that 
the fittings of the numerical results for P show that the 
assymptotic behaviors of this magnitude for random ma- 
trices and for the present billiard are the same, and that 
significant differences between the two models are only 
found for small L. These results would suggest that, for 
L — > oo, whereas averaged properties of quantum chaotic 
billiards approach those of random matrices, fluctuations 
are much larger in the former. 
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FIG. 4. Scaling behavior of the relative fluctuation of the 
participation ratio for random matrices of the GOE (squares) 
and for our quantum billiard (circles). Fits are shown as 
dashed lines. 

The results for P allow us to connect with a question of 
much recent interest. We refer to the eventual multifrac- 
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tality jl8],[l9) of the wavefunctions predicted |^0|,^l| and 
numerically calculated [ fl6| for prelocalized states in dis- 
ordered two dimensional systems. The point is whether 
this exotic behavior of the eigenstates could also be a 
characteristic of chaotic wavefunctions. Multifractality 
occurs whenever t qa oc L~ T ( q \ r(q) being a non-integer; 
in particular r{q) = (q — l)D(q), where the D(q) are the 
generalised fractal dimensions. Our results for P a (or 
t 2a ) averaged over the whole energy range, indicate that 
this behavior is not expected in the present from 
the above fittings it is concluded that D(2) = 2. In or- 
der to ensure this conclusion we have carried out several 
analysis in selected energy ranges by means of the above 
method and by the standard box-counting method fl~9| l 
for clusters of a fixed size. If size effects are properly ac- 
counted for, all results point to the same conclusion: the 
wavefunctions of the billiard herewith investigated do not 
show multifractal behavior. In fact all results for t qa L 2q 
can be most accurately fitted by parabollic functions. 

Summarizing, we have presented a new model of quan- 
tum chaotic billiards which is an efficient implementation 
of surface roughness. The essential feature of the model 
is the inclusion of diagonal disorder at the surface of the 
system. The spectral statistics of this billiard changes 
through the band in a manner not previously reported 
in other models of chaotic billiards. In particular, ex- 
ponentially localized, quasi-ideal, surface resonances and 
chaotic states are found to exist within the band. The 
probability amplitude of chaotic eigenstates is reminis- 
cent of that found in more standard chaotic billiards and 
for prelocalized states in two dimension. We have also 
shown that whereas the assymptotic behavior of the the 
participation ratio in the present billiard is almost iden- 
tical to that found in random matrices, the standard de- 
viation (fluctuations) of that magnitude is much larger in 
the former. This results suggests that mapping between 
chaotic billiards and the random matrix problem should 
only be expected for averaged properties and not for their 
fluctuations. It is very likely that these features are not 
exclusive of the present billiard, and that the behavior 
of quantum chaotic billiards cannot be fully described, 
at least at a mesoscopic level, by random matrices. Fi- 
nally, we note that the simplicity of our model allows the 
study of several situations of physical interest including 
the case of 3D billiards. 
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